Could pressureless dark matter have pressure? 
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Abstract 

A two-fluid dark matter model, in which dark matter is represented as a two-component fluid 
thermodynamic system, without interaction between the constituent particles of different species. 
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^J^l and with each distinct component having a different four- velocity, was recently proposed in Harko & 

Lobo, [Phys. Rev. D83, 124051 (2011)]. In the present paper we further investigate the two-fluid 
. dark matter model, by assuming that the two dark matter components are pressureless, non- 

I comoving fluids. For this particular choice of the equations of state the dark matter distribution 

' can be described as a single anisotropic fluid, with vanishing tangential pressure, and non-zero 

■ radial pressure. We investigate the properties of this model in the region of constant velocity 

galactic rotation curves, where the dynamics of the test particles is essentially determined by the 
- ^ _ dark matter only. By solving the general relativistic equations of mass continuity and hydrostatic 

H ■ equilibrium we obtain the geometric and physical parameters of the dark matter halos in the 

constant velocity region in an exact analytical form. The general, radial coordinate dependent, 
functional relationship between the energy density and the radial pressure is also determined, and 
it differs from a simple barotropic equation of state. 

PACS numbers: 04.50.Kd, 04.20. Cv, 04.20.Fy 



* Electronic address: lharko@hkucc.hku.hkl 
^Electronic address: | flobo@cii.fc.ul.pt" 



1 



I. INTRODUCTION 



The Concordance Cosmological Model, usually referred to as the A cold dark matter 
(ACDM) model, has proved to be very successful in explaining cosmological observations 
across a wide rage of length scales, from the cosmic microwave background (CMB) anisotropy 
to the Lyman-a forest [ll, |2|. In this model, nonbaryonic collisionless cold dark matter makes 
up to 23% of the total mass content of the Universe. In the ACDM model, dark matter 
consists of cold neutral weakly interacting massive particles, beyond those existing in the 
Standard Model of Particle Physics. However, up to now no dark matter candidates have 
been detected in particle accelerators or in direct and indirect searches. Many particles 
have been proposed as possible candidates for dark matter, the most popular ones being the 
Weakly Interacting Massive Particles (WIMP) and the axions (for a review of the particle 
physics aspects of dark matter see j^). The interaction cross section of dark matter particles 
with normal baryonic matter is assumed to be extremely small. However, it is expected to 
be non-zero, and therefore the direct experimental detection of dark matter particles may 
be possible. Superheavy particles, with mass > 10^° GeV, have also been proposed as dark 
matter candidates. But in this case observational results show that these particles must 

r 

either interact weakly with normal matter, or they must have masses above 10^^ GeV j^. 
Scalar field models, or other long range coherent fields coupled to gravity have also been 
proposed to model galactic dark matter [5]-[2l|. The possibility that dark matter could 
be described by a fluid with non-zero effective pressure was also investigated 22[ 23 1. In 
particular, it was assumed in 2J| that the equation of state of the dark matter halos is 
polytropic. The fit with a polytropic dark halo improves the velocity dispersion profiles. 
The possibility that the galactic dynamics of massive test particles may be understood 
without the need for dark matter was explored in the context of modified theories of gravity 



m 



25j-|32|. 



On galactic scales observational data seem to disagree with the ACDM model predictions. 
High resolution N-body simulations have shown that the predicted number of subhalos is 
an order of magnitude larger than what has been observed 33|]. Another discrepancy arises 
when comparing the density profiles of dark halos predicted in simulations with those derived 
from observations of dwarf spheroidal (dSph) galaxies and Low Surface Brightness galaxies 
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(LSBs). N-body simulations predict an universal cuspy density profile |34j. |35| 



PNFw{.r) ^ ^\ — — (1) 

where is a scale radius and ps is a characteristic density. On the other hand observations 
based on high-resolution rotation curves show, instead, that the actual distribution of dark 
matter is much shallower than the above, thus indicating that a cored halo is preferred in 
an important fraction of low-mass galaxies jsGj, 

PB^r) ^ y^^^l^^y (2) 

where tq is the core radius and po is the central density. The observational Burkert density 
profile, given by Eq. ([2]), resembles an isothermal profile in the inner regions, i.e., r ^ ro, 



predicts a finite central density, po? and leads to a mass profile that diverges logarithmical 



for increasing r, which is consistent with cosmological cold dark matter predictions [3j, |35|. 

These discrepancies between theoretical predictions and observations might be overcome 
by considering other alternative dark matter models. The possibility that dark matter is 
a mixture of two non-interacting perfect fluids, with different four-velocities and thermo- 



dynamic parameters, was proposed recently in 371]. By introducing a rotation of the four- 



velocity vectors the two-fluid model can be reduced to an effective single anisotropic fluid 



model, with distinct radial and tangential pressures |38H40l| . By assuming a non-relativistic 
kinetic model for the dark matter particles, the density profile and the tangential velocity 
of the dark matter mixture have been obtained by numerically integrating the gravitational 
field equations. The cosmological implications of the model have also been briefly ana- 
lyzed, and it was shown that the anisotropic two-fluid model isotropizes in the large time 
limit. Two fluid dust models have also been considered in a general relativistic framework 



m 



41 



It is the purpose of the present paper to further investigate the idea proposed in |37| . 
by considering the specific case of two, non-interacting, pressureless dark matter fluids. For 
this configuration the model reduces to a single anisotropic fluid, with vanishing tangential 
pressure. We investigate the properties of this model in the region of constant galactic 
velocity rotation curves, where the solution of the basic equations can be obtained in an 
exact analytical form. Thus, the radial coordinate dependence of all relevant geometric and 
physical parameters of the dark matter halos is explicitly determined. In particular, we 



obtain the general, r- dependent, functional relationship between the energy density and the 



rom the simple barotropic equation of state 
M. 
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radial pressure of the dark matter, which differs 
previously considered in the physical literature 

The present paper is organized as follows. The two-fiuid model of the dark matter halos 
is briefly reviewed in Section HIl The general relativistic structure equations for anisotropic 
fluids are written down in Section llllt and the tangential velocity of test particles in stable 
circular orbits is obtained as a function of the geometric metric tensor. The general solution 
of the gravitational field equations in the constant velocity region of the dark matter halos 
is obtained in Section HVl We discuss and conclude our results in Section |Vl In the present 
paper we use the natural system of units with c = G = h = 1. 



II. DARK MATTER AS A MIXTURE OF TWO PERFECT FLUIDS 

We start our study of dark matter by assuming that it consists of a mixture of two perfect 
fluids, with energy densities and pressures pi, pi and p2, P2, respectively, and with four 
velocities f/^ and VV^, respectively. The fluid is described by the total energy- momentum 
tensor T^"^ , given by 

T"" = (pi + pi) U'-U" - Pig''" + (p2 + P2) WW - p^g^". (3) 

The four- velocities are normalized according to 1/^11^ = 1 and W^W^ = 1, respectively. 
The study of the physical systems described by an energy-momentum tensor having the 
form given by Eq. can be significantly simplified if we cast it into the standard form of 
perfect anisotropic fluids. This can be done by means of the transformations If^ — )• U*'^ and 



W*^', respectively, so that |38|-|40| 




cos a ^ P2+V2 gjj^ Q, 
pi+pi 



^^4^ sin a cos a 

P2+P2 




(4) 



representing a "rotation" of the velocity four-vectors in the {U^, W^) velocity space. No- 
tice that the transformations given by Eqs. (jll) leave the quadratic form (pi + pi) U^V^ + 
(P2 + P2) W^W^ invariant. Thus we have 

T"'' {U,W) = T"" {U\W*) . (5) 
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As for the vectors U*^ and W*^ we assume that one is timehke, while the other is spacehke, 



so that U*W, = 0. 



With the use of the latter relationship and Eqs. ([!]) we obtain the rotation angle as 

Pi + Pi - {P2+P2) 



(6) 



Next we define the following quantities 38|-|40[| : 



(7) 



e = T^-^V.V, = (pi + pi) f/*"f/: - (pi + P2) , (8) 
a = T^'x.Xu = (pi + P2) - (P2 + P2) W*^W:, (9) 

n = pi+P2, (10) 

respectively. Thus, the energy-momentum tensor of the two non-interacting perfect fiuids 
can be written as 

T^" = {e + Tl) V^V - n^^" + (a - n) x'^x", (n) 

where 

V^V^ = 1 = -x^Xm, (12) 

and 

X^V; = 0. (13) 

Note that the energy-momentum tensor given by Eq. ffTTl) is the standard form for 
anisotropic fiuids [4^. 

The energy density e and the radial pressure a are given by 



£ = -{pi+p2-pi -P2) + 



(pi + Pi + P2 + P2f + 4 (pi + pi) (p2 + P2) [{U^'W^,f - 1] ^ 



(14) 



and 



^ = ~2 (^1 + ^2 - Pi - P2) + 



^\/(Pi + Pi - P2 - P2)' + 4 (pi + pi) (p2 + P2) {Um^f, 



(15) 



respectively 



38|-|4Q|. 



In comoving spherical coordinates 
= V^ = V'^ = 0, V^Vo = 1, and x° - 



and = 6 we may choose 



X 



0, x^Xi = -1 |38|-|4Q|. Therefore the 



components of the energy-momentum of two non-interacting perfect fluids take the form 



-^0 —=^5-^1 



Ti = -n, 



(16) 



where e is the total energy-density of the mixture of fluids, a = Pr is the pressure along the 
radial direction, while IT = Pi is the tangential pressure on the r = constant surface. 



III. ANISOTROPIC FLUIDS IN SPHERICALLY SYMMETRIC STATIC SPACE- 
TIMES 

In the following we restrict our study of the two-component dark matter to the static 
and spherically symmetric case, with the metric represented as 

ds^ = e'^'^dt^ - e^^'^dr^ - [dd^ + sin^ M<P^) . (17) 

For the metric given by Eq. f|T7|) . the Einstein gravitational field equations, describing 
the mixture of two pressureless fluids, take the form 

3 - - ) + 3 = 8^^' (18) 
?7r(T, (19) 

and 

a' 4U-a 
v' = -2 + , 21 

e + a r e + a 

where ' = d/dr. Equation ( l2Ti) is the consequence of the conservation of the energy- 
momentum tensor, T^^,;^ = 0. Eq. ( JTSl) can be easily integrated to give 

e-^ = l-'-Mn, (22) 

r 

where M(r) = 4tt J er'^dr. With the use of Eqs. (fT9|) and fl2T|) we obtain the mass continuity 
equation, and the generalized Tolman-Oppenheimer-Volkoff (TOV) equation describing the 
static pressureless two-fluid dark matter configurations, 

dM 
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A' 




r 







dr 
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ATier\ (23) 



da ^ je + a) {ATrar^ + M) , 2 . 

dr r2(l-2M/r) ^ ^ 

The galactic rotation curves provide the most direct method of analyzing the gravitational 
field inside a spiral galaxy. The rotation curves are obtained by measuring the frequency 
shifts z of the 21-cm radiation emission from the neutral hydrogen gas clouds. Usually 



astronomers report the resulting z in terms of a velocity field Vtg |43l-l46| . The line element 
given by Eq. (1171) can be rewritten in terms of the spatial components of the velocity, 
normalized with the speed of light, measured by an inertial observer far from the source, as 
ds'^ = (it^ (1 — f^), where 



I dr\^ 2 f dQ \ ^ 



(25) 



For a stable circular orbit dr/dt = 0, the tangential velocity of the test particle can be 
expressed as 



<^e-v(f ) , (26) 



where dQ"^ = dd'^ + sin^ dd^)^. 

In terms of the conserved angular momentum / = and energy E = e'^^'^H, the tan gen- 
tial velocity is given, for = 7r/2, by 

< = "^w 

By taking into account the explicit expressions for / and E we obtain for the tangential 
velocity of a test particle in a stable circular orbit the expression 

< = Iru'. (28) 

The tangential velocity profile can be expressed as a function of the total density and of 
the mass of the dark matter from Eq. (fT9|) as 



r, , , 47r o-r^ + M 
= r(l-2M/ry 



The TOV Eq. (!24|) . can be written in the following equivalent form 

da (e + a) o , , 2 ,^ , , . 
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IV. DARK MATTER AS A MIXTURE OF TWO PRESSURELESS FLUIDS 



In order to describe dark matter we adopt the kinetic model considered in [37|]. The 
energy-momentum tensor Tj^ associated to the frozen distribution of dark matter is given by 
'^u ~ 9 1 d^Pf{p)p'^Pu/p^, where f{p) is the dark matter particle distribution function, is 
the four-momentum, p is the three-momentum, with absolute value p, and g is the number 
of helicity states, respectively. The energy density e of the system is defined as 

e = 1 1 Ef{p)d'p, (31) 

while the pressure of a system with an isotropic distribution of momenta is given by 

P = 1 1 pvfip)d'p ^fip)d'p, (32) 



where the velocity v is related to the momentum hj v = p/E 47|, |48|]. In the non-relativistic 



regime, when E ^ m and p ^ mv, the density p of the dark matter is given by p = mn 



47 



48| 



where n is the particle number density, while its pressure P can be obtained as 

P = U T^f(P)d'P ^ 47r| / P-dp, (33) 
3 J E 3 J m 

giving 

P=^P, (34) 

where {v'^) is the average squared velocity of the particle, and (?T^)/3 is the one- dimensional 
velocity dispersion. In the non-relativistic approximation given by Eqs. ( 133|) and the 
velocity dispersion is a constant only for the case of the non-degenerate ideal Maxwell- 
Boltzmann gas whose Newtonian analogue is the isothermal sphere. In the following we 
will consider only the specific case of the ideal gases in the non-relativistic regime. Since 
for non-relativistic particles (w^)/3 <^ 1, and their energy density is also very small, we 
can consider, with a very good approximation, that each component of the dark matter is 
pressureless, P ^ 0. Therefore we assume that dark matter consists of the mixture of two 
pressureless fluids, so that Pi = P2 ~ 0. The two fluids have different four-velocities, and 
hence we may write 

W^W, = 1 + ^, (35) 

where in the general case b is an arbitrary function of the radial coordinate r. The functional 
form of b can be estimated by using Eq. ([6]). If the angle a between the two four- velocities 
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is small, we obtain 

W^W, ^ 4^"- (36) 
VP2/P1 

In the following we assume that the number density of the two components is roughly the 
same, rii = n2 = n, but the dark matter particles have different masses mi 7^ 7712- Thus 
P2/P1 = /mi, and b can be estimated as 

1 -ms/mi \ 

/ , " - 1 . (37) 
A/ma/mi J 

The components of the effective anisotropic energy-momentum tensor for this mixture of 
fluids is given by 



^ = 2 (pi + P2) + 2 V (-"1 + p-^y + + p^p^^ (3s) 

^ = = (Pl + P2) + ^7(Pl-p2)' + (& + 2)VlP2, (39) 

n = = 0, (40) 

respectively. Therefore the mixture of the two pressureless fluids can be modeled as an 
anisotropic fluid with vanishing tangential pressure. In the following we will restrict our 
study to the constant tangential velocity region, so that we can take Vtg = constant. This 
assumption determines, with the use of Eq. ( l28i) . the functional form of the metric tensor 
component in the dark matter region as 

where rQ^y is an arbitrary integration constant. The gravitational field equation Eq. (!20|) can 
be written as 



A' = (42) 



which immediately gives 



e 



where tqa is an arbitrary integration constant. The density profile of the two component 
dark matter mixture follows from Eq. (fTSl) . and is given by 



1 



1 _ *^ 



1 + ^tg V'^OA 



(43) 
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while the effective pressure of the two component dark matter is given by 



Sua 



-2vl/{l+vl) 



- 1 



(44) 



Equation (jHj) enables us to introduce the radius R of the dark matter distribution, defined 
as the vacuum boundary of the system, where the radial pressure vanishes, a {R) = 0. This 
condition gives 

ro.= (l + 2<)'(^^^-)/^^-i?. (45) 
Thus, the mass profile of the dark matter is provided by 



m(r) 



r - rox 



rox 



:i-?.)/(i+-?«) 



The total mass of the dark halo is given by M = m{R), and, by using the relation 



\roxJ 

can be obtained in a simple form as 



R \ 



l + 2vl \roxJ 



M 



The dark matter density at the boundary surface is given by 



Svre {R) = 



(46) 



(47) 



(48) 



(49) 



The equation of state of dark matter can be formally obtained from Eq. (l30l) . Since 
n = 0, Eq. ( !30l) can be written as a first order differential equation. 



da 
dr 



2 + < 



which can be integrated to give 

a(r) = r-(2+"*«) 



o"o + / r^^^^ae (r) dr 



(50) 



(51) 



where cxo is an arbitrary integration constant. 

However, in the present model the equation of state is not of a simple barotropic form, and 
the functional relation between e and a is also r-dependent. Once e and a are known, from 
Eqs. ffl^ and (ITSl) one can obtain the density profiles of the two dark matter components, 
which are given by 



Pi 



2 (8 + 6 {h + A))ea 



a" 



1,2. 



(52) 
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V. DISCUSSIONS AND FINAL REMARKS 



In the present paper we have further considered the theoretical possibihty, proposed in 



371], that dark matter may be modeled as a mixture of two non-interacting perfect fluids, 



with different four-velocities. The most interesting characteristic of this model is that it is 
formally equivalent to a single anisotropic fluid. It is important to note that this equivalence 
is general, and it is independent of the stationarity condition of the metric assumed in the 
present paper. Mathematically, this result follows from the invariance of the total energy- 
momentum tensor of a two-fluid system with respect to the group of rotations in the velocity 
space. Thus, the same approach can be used for the study of time-dependent dark matter 
halos, like, for example, those resulting from the collision of two galaxies. In particular, for 
the specific case of two non-interacting pressureless fluids, the two fluid model reduces to a 
single anisotropic fluid with vanishing tangential pressure. 

By considering the flat galactic rotation curves region, we have completely solved the 
general relativistic structure equations of the pressureless two fluid dark matter halos. We 
have found a general, r-dependent functional relationship between the energy density and the 
radial pressure of the dark matter halos, which differs radically from the simple barotropic 
equation of state previously considered in the literature. On the other hand, in was shown 



m 



24| that half of the 14 galaxies considered in the study are well fit by the polytropic halo 



model, despite its serious physical simplifications. 

The necessity of considering non-standard dark matter models with pressure is justified 
by our uncertainties in the knowledge about the nature of the dark matter particles, as well 
as by the fact that this model gives a much better description of the observational results, 
as compared to the pressureless case. On the other hand, the very successful standard 
ACDM cosmological model is based on the fundamental assumption of the existence of the 
pressureless dark matter. 

The apparent contradiction between these two approaches could be solved by assuming 
that dark matter is a mixture of two pressureless fluids, with different four-velocities. Then 
dark matter can be described as a single, anisotropic fluid. Similarly to the case of the 
NFW profile, the anisotropic dark matter profile given by Eq. ( H3|) diverges for r — )■ 0, thus 
predicting a cuspy profile. However, this unphysical divergence of the density is the result 
of our assumption of a constant tangential velocity, which enables us to obtain a complete 
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analytical solution of the full set of the gravitational field equations. By assuming an ar 



bitrary velocity profile, the field equations can be solved numerically, as done in [37|, and 
the obtained density distributions have a non-singular, core-like profile. From an analytical 
point of view, one way to deal with this unphysical feature would be to perform a smooth 
matching between the anisotropic fluid distribution, and a small central region with a regular 
density profile. This approach would lead to a better analytical description of the dark mat- 
ter halos. Another possibility would be to assume that indeed each of the two components 
of the dark matter fluid have some pressure 37|, and consider some approximate solutions 
of the field equations. Some of these possibilities, as well as their impact on the dark matter 
distribution and properties, together with the cosmological implications, will be considered 
in a future study. 
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